Higher-order gravity and the cosmological background of gravitational waves 
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The cosmological background of gravitational waves can be tuned by the higher-order corrections 
to the gravitational Lagrangian. In particular, it can be shown that assuming i? 1+e , where e indicates 
a generic (eventually small) correction to the Hilbert-Einstein action in the Ricci scalar R, gives a 
parametric approach to control the evolution and the production mechanism of gravitational waves 
r***» . in the early Universe. 
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Several issues coming from Cosmology and Quantum Field Theory suggest to extend the Einstein General Relativity 
in order to cure shortcomings emerging from observations and self-consistent unification theories. In early time 
Cosmology, the presence of Big Bang singularity, flatness and horizon problems [l[ led to the result that Standard 
Cosmological Model @, is inadequate to describe the Universe at extreme regimes. On the other hand, General 
Relativity is a classical theory which does not work as a fundamental theory, when one wants to achieve a full 
I ' quantum description of spacetime (and then of gravity). Due to these facts and, first of all, to the lack of a self- 
consistent Quantum Gravity theory, alternative theories of gravity have been pursued in order to attempt, at least, 
a semi-classical scheme where General Relativity and its positive results could be recovered. A fruitful approach 
has been that of Extended Theories of Gravity which have become a sort of paradigm in the study of gravitational 
interaction based on corrections and enlargements of the Einstein scheme [!, |j| . 

Besides fundamental physics motivations, these theories have acquired a huge interest in cosmology due to the fact 
that they "naturally" exhibit inflationary behaviors able to overcome the shortcomings of Standard Cosmological 
Model (based on General Relativity). The related cosmological models seem very realistic and, several times, capable 
. of matching with the observations jf| [||. Recently, Extended Theories of Gravity are going to play an interesting 
role to describe the today observed Universe. In fact, the amount of good quality data of last decade has made it 
possible to shed new light on the effective picture of the Universe. Type la Supernovae (SNela) @, anisotropics 
£^ in the cosmic microwave background radiation (CMBR) [8], and matter power spectrum inferred from large galaxy 
surveys [§] represent the strongest evidences for a radical revision of the Cosmological Standard Model also at recent 
!* ■ epochs. In particular, the concordance ACDM model predicts that baryons contribute only for ~ 4% of the total 
. £h ! matter - energy budget, while the exotic cold dark matter (CDM) represents the bulk of the matter content (~ 25%) 
^ . and the cosmological constant A plays the role of the so called "dark energy" (~ 70%) 0. Although being the best 
fit to a wide range of data [ll[, the ACDM model is severely affected by strong theoretical shortcomings [12] that 
have motivated the search for alternative models (l3| . Dark energy models mainly rely on the implicit assumption 
that Einstein's General Relativity is the correct theory of gravity indeed. Nevertheless, its validity on the larger 
astrophysical and cosmological scales has never been tested and it is therefore conceivable that both cosmic 
speed up and dark matter represent signals of a breakdown in our understanding of gravitation law so that one should 
consider the possibility that the Hilbert-Einstein Lagrangian, linear in the Ricci scalar R, should be generalized 
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From a genuine astrophysical viewpoint, Extended Theories of Gravity do not urgently require to find out candidates 
for dark energy and dark matter at fundamental level (till now they have not been found!); the approach starts from 
taking into account only the "observed" ingredients (i.e. gravity, radiation and baryonic matter); it is in full agreement 
with the early spirit of General Relativity which could not act in the same way at all scales (for comprehensive reviews 
on the argument, see [20L [2l|). In fact, General Relativity has been successfully probed in the weak filed limit (e.g. 
Solar System experiments) and also in this case there is room for alternative theories of gravity which are not at all 
ruled out, as discussed in several recent studies [H, [H, [H, HE] . In particular, it is possible to show that several f(R) 
models could satisfy both cosmological and Solar System tests [26[ , could be constrained as the scalar-tensor theories 
(22l [27] | , could give rise to new effects capable of explaining the so called Pioneer anomaly (see for example [28| and 
the references therein). 

However, also considering the recent interest in the argument, a comprehensive effective theory of gravity, acting 
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consistently at any scale, is far, up to now, to be found out, and this demands an improvement of observational datasets 
and the search for experimentally testable theories. A more pragmatic point of view could be to "reconstruct" the 
suitable theory of gravity starting from data. The main issues of this "inverse " approach is matching consistently 
observations at different scales and taking into account wide classes of gravitational theories where "ad hoc" hypotheses 
are avoided. In principle, the most popular dark energy models can be achieved by considering f{R) theories of gravity 
[29l I30I ] and the same track can be followed, at completely different scales, to match galactic dynamics [3l[. This 
philosophy can be taken into account also for the cosmological stochastic background of gravitational waves (GW) 
which, together whit cosmic microwave background radiation (CMBR) [32j ]. would carry, if detected, a huge amount 
of information on the early stages of the Universe evolution. However, also in this case, a key role for the production 
and the detection of this relic graviton background is played by the adopted theory of gravity which gives rise to 
specific early Universe models [H, H3, HU . 

In this letter we want to face the problem of how a theory of gravity with a Lagrangian of the form R 1+e , where 
e is a deviation (e.g. small) with respect to General Relativity, could be related to the cosmological background of 
GWs. They are perturbations /i M „ of the metric g^ u which transform as 3-tensors. Following [2j, the GW-equations 
in the transverse-traceless gauge are 

Uh\ = (1) 

where □ = {—g)^ 1 ^ 2 d ll (~g) 1 ^ 2 g^ v di, is the d'Alembert operator; these equations are deduced from the weak-field 
limit, in vacuum, of the Einstein field equations. The Latin indexes run from 1 to 3, the Greek ones from to 3. Our 
task is now to derive the analog of Eqs. |T]) assuming a correction to the Hilbert-Einstein action given by 

•A = I <PxV=gf R 1+e ■ (2) 
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From now on we will assume units for which fo = 1. It is easy to show that 

i? 1+£ = R-R e ~i?(l + (logi?)e + 0(e 2 )) ~ R + eR log R, (3) 

and this form of the action can be assumed: i) to define the right physical dimensions of the coupling constant; ii) 
to control the magnitude of the corrections with respect to the standard Einstein gravity. Assuming a conformal 
transformation, the extra degrees of freedom related to the higher order gravity can be recast into a physically 
significant scalar field (for a detailed discussion see [H, H3] and references therein) 

SV = e 2< V with e 2 * = ^ = /'(JJ)~(l + e + elogJ2) (4) 

where the prime indicates the derivative with respect to the Ricci scalar R and $ is the "conformal scalar field" . The 
conformally equivalent Hilbert-Einstein action is 

* = h I ^ d ' x [k +£ ( $,$ ^)] (5) 

where C ($,$;^J is the conformal scalar field Lagrangian derived from 

R = e- 2 * (R - 6D$ - 6$ ;5 $ ;<5 ) . (6) 

Deriving the field equations from ([5]), the GW-equations are 

Dhl = (7) 

expressed in the conformal metric <7 M „. Since no scalar perturbation couples to the tensor part of gravitational waves, 
we have 

H = g ll 5g a = e- 29 g l h^5 9il = h\ (8) 

which means that h\ is a conformal invariant. 

As a consequence, the plane-wave amplitudes h^t) = h(t)e? exp(ik m x m ), where is the polarization tensor, are 
the same in both metrics. This fact will assume a key role in the following discussion. The d'Alembert operator 
transforms as 

□ = e~ 2 * (□ + 2$ ;A 9 ;A ) (9) 
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and this means that the background is changing with frame while the tensor plane- wave amplitudes not. 

In order to study the gravitational waves in the cosmological background, the operator @ has to be specified for 
a Friedmann-Robertson- Walker metric and then Eq. ([7]) becomes 

h+ (W + 2$) h + k 2 a' 2 h = (10) 

being a(t) the scale factor, k the wave number and h depending only on the surviving component of h\ in the 
homogeneous and isotropic metric. 

It is worth stressing that Eq. (|10p applies to any f(R) theory whose conformal transformation can be defined as 
e 2 * = f'(R)- The GW amplitude depends on the specific cosmological background a(t) and the specific theory of 
gravity which can be parameterized by e in our case ©. For example, assuming time power law behaviors for a(t) 
and f'(R(t)), we have 

f'(R) = f (t/t Q ) m , a(t) = a (t/t ) n . (11) 
vn 

For a generic f(R) = R 1+e , it is e = — — . Eq. (TIT)]) can be recast in the form 

h + (3ra - 2e) r v h + k 2 a {t /t)) 2n h = (12) 

whose general solutions are J Q 's ar Bessel functions. In Fig. ([I]) some examples are given. The plots are labelled by 
the set of parameters {m, n, e} which assign the time evolution of &(t) and a(t) with respect to a given power-law 
theory. 

The time units are in terms of the Hubble radius H" 1 . It is clear that the conformally invariant plane- wave 
amplitude evolution of the tensor GW strictly depends on the cosmological background "tuned" by e. 

Let us now take into account the issue of the production of GWs contributing to the cosmological stochastic 
background. Several mechanisms can be considered as vacuum fluctuations, phase transitions [33J, cosmological 
populations of astrophysical sources (38j and so on. In principle, we could seek for contributions due to every high- 
energy process in the early phases of the Universe evolution. 

It is important to distinguish processes coming from transitions like inflation, where the Hubble flow emerges in 
the radiation dominated phase and processes, like the early star formation rates, where the production takes place 
during the dust dominated era. In the first case, stochastic GW background is strictly related to the cosmological 
model. This is the case we are considering here which is strictly connected to the specific theory of gravity. In 
particular, one can assume that the main contribution to the stochastic background comes from the amplification of 
vacuum fluctuations at the transition between an inflationary phase and the radiation dominated era. However, in 
any inflationary model, we can assume that the GWs generated as zero-point fluctuations during the inflation undergo 
adiabatically damped oscillations (<~ 1/a) until they reach the Hubble radius H~ l . This is the particle horizon for the 
growth of perturbations. On the other hand, any previous fluctuation is smoothed away by the inflationary expansion. 
The GWs freeze out for a/k 3> H^ 1 and reenter the radius after the reheating in the Friedmann era (see also 
[3{||4(|). The reenter in the radiation-dominated or in the dust-dominated era depends on the scale of the GW. After 
the reenter, GWs can be detected by their Sachs- Wolfe effect on the temperature anisotropy AT/T at the decoupling 
[4lj . If <E> acts as the inflaton we have $ <C H during the inflation. Adopting the conformal time dr\ = dt/a, Eq. (jTUJ) 
reads 

h" + 2— ti + k 2 h = (13) 

X 

where \ = ae*. The derivation is now with respect to rj. An inflationary behavior is achieved for a(t) = aoexp(Ht), 
then i] — J dt/a = (aH)^ 1 and x'/x = ~ V 1 - The exact solution of (fl"3"|) is 

h(rj) = k~ z l 2 \f2jk [Ci (sinfery - cosfc??) + C 2 (sinfcr; + cosfcr;)] (14) 

Inside the radius H -1 , we have krj ^> 1. Considering the absence of gravitons in the initial vacuum state, we have 
only negative-frequency modes and then the adiabatic behavior is 

h = k l ' 2 J2hi— C&M-ikrf). (15) 
an 

where C is the parameter to be tuned later. At the first horizon crossing (aH — k), the averaged amplitude Ah = 
(k/2n) 3/2 \h\ of the perturbatio n is 

A h = ^C. (16) 
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Figure 1: Qualitative evolution of the GW amplitude for some choices of a(t), 3>(i) and R +E . Time-scales, amplitudes and 
"chirps" of GWs strictly depend on the value of e as it is especially shown in the last two (unrealistic) cases (e extremely large). 

When the scale a/k becomes larger than the Hubble radius H , the growing mode of evolution is constant, i.e. it 
is frozen. This situation corresponds to the limit — kr\ <C 1 in Eq. ([HI) . Since <E> acts as the inflaton field, we have 
$ ~ at reenter (after the end of inflation) . Then the amplitude Ah of the wave is preserved until the second horizon 
crossing; after it can be observed, in principle, as an anisotropy perturbation on the CMBR. It can be shown that 
AT/T < Ah, as an upper limit to Ah, since other effects can contribute to the background anisotropy. From this 
consideration, it is clear that the only relevant quantity is the initial amplitude C in Eq. (|15p , which is conserved until 
the reenter. Such an amplitude depends on the fundamental mechanism generating perturbations. Inflation gives rise 
to processes capable of producing perturbations as zero-point energy fluctuations. Such a mechanism depends on the 
gravitational interaction and then (AT/T) could constitute a further constraint to select a suitable theory of gravity. 
Considering a single graviton in the form of a monochromatic wave, its zero-point amplitude is derived through the 
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commutation relations: 



[h{t,x), Tv h (t,y)] = iS 3 (x-y) 



(17) 



calculated at a fixed time t, where the amplitude h is the field and 717, is the conjugate momentum operator. Writing 
the Lagrangian for h 



in the conformal FRW metric (recall that the amplitude h is conformally invariant), we obtain 

dC 



ir h = ^ = e 2 *a 3 h 

dh 



Eq. (fT7|) becomes 



h{t,x), h{y,y) 



.d 3 (x-y) 



and the fields h and h can be expanded in terms of creation and annihilation operators 

1 



h(t, x) 



(27r) 3 / 2 



d 3 k [h(t)e~ lkx + h*(t)e 



-\-ikx~\ 



(18) 



(19) 



(20) 



(21) 



h(t, x) 



1 



d 3 k 



(27r) 3 / 2 

The commutation relations in conformal time arc then 



h{t)e~ lkx +h*(t)e 



-\-ikx 



(22) 



(23) 



Using (|15p and (HHJ), we obtain C — \/2TT 2 He *, where H and <E> are calculated at the first horizon-crossing and 



or, in general, 



2 2 V 2 



^i?(l-|logi?+^(logi?) 2 



H 



(24) 



(25) 



Clearly the amplitude of GWs produced during inflation depends on the given theory of gravity that, if different 
from General Relativity, gives extra degrees of freedom capable of affecting the cosmological dynamics. On the other 
hand, the Sachs- Wolfe effect related to the CMBR temperature anisotropy could constitute a test for the theory of 
gravity at early epochs, i.e. at very high redshift. This probe could give further constraints on the GW-stochastic 
background, if Extended Theories of Gravity are independently probed at other scales. In summary, assuming that 
primordial vacuum fluctuations produce GWs, beside scalar perturbations, kinematical distortions and so on, the 
initial amplitude of these ones is a function of the assumed theory of gravity and then the stochastic background can, 
in a certain sense, be "tuned" by the theory. Conversely, data coming from the Sachs- Wolfe effect could contribute 
to select a suitable e which can be consistently matched with other observations. If eventually, one finds e rigorously 
zero, this could be a further test for General Relativity also at very early cosmological scales. This result could be 
obtained very soon through space and ground based interferometers (42|. 
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